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We solve the problem of exact minimization of the Lawrence-Doniach (LD) free-energy functional 
in parallel magnetic fields. We consider both the infinite in the layering direction case (the infinite 
LD model) and the finite one (the finite LD model). We prove that, contrary to a prevailing view, 
the infinite LD model does not admit solutions in the form of isolated Josephson vortices. For 
the infinite LD model, we derive a closed, self-consistent system of mean-field equations involving 
only two variables. Exact solutions to these equations prove simultaneous penetration of Josephson 
vortices into all the barriers, accompanied by oscillations and jumps of the magnetization, and yield 
a completely new expression for the lower critical field. Moreover, the obtained equations allow us 
to make self-consistent refinements on such well-known results as the Meissner state, Fraunhofer 
oscillations of the critical Josephson current, the upper critical field, and the vortex solution of 
Theodorakis [S. Theodorakis, Phys. Rev. B 42, 10172 (1990)]. Our consideration of the finite 
LD model illuminates the role of the boundary effect. In contrast to the infinite case, an explicit 
analytical solution to the Maxwell equations of the finite case does not preclude the existence of 
localized Josephson vortex configurations. By the use of this solution, we obtain a self-consistent 
description of the Meissner state. Finally, we discuss some theoretical and experimental implications. 

PACS numbers: 74.80.Dm, 74.20.De, 74.50. +r 



I. INTRODUCTION 

We obtain exact analytical solutions to the phenomenological Lawrence-Doniach (LD) model for layered super- 
conductors in external parallel magnetic fields. We consider both the infinite in the layering direction case (the infinite 
LD model) and the finite one (the finite LD model). This paper should be considered as a logical continuation of our 
previous study of layered superconductors on the basis of a microscopic approach . 

At present, the LD model is widely used for the description of low-T c layered superconductors and superlattices as 
well as high-T c superconductors exhibiting the intrinsic Josephson effect |3|,|4|] . Surprisingly, despite a large number of 
theoretical publications on this subject, it has not been realized yet that the problem of the parallel magnetic field is 
exactly solvable. Up to now, actual analytical solutions with different degrees of accuracy have been obtained only 
for relatively simple particular cases of the infinite LD model: the Meissner state ||, Fraunhofer oscillations of the 
critical Josephson current B, the upper critical field H C 2(T) and the vortex state Q in the intermediate field 

regime. 

Unfortunately, the calculations of the lower critical field H c i ||[l0| , based on the assumption of isolated Josephson 
vortex penetration, raise questions. In these calculations, one employs an anisotropic continuum approximation outside 
the so-called "Josephson vortex core region" ]TT| , completely neglecting the intrinsic discreteness of the LD model. 
As has been recently shown by Farid JL2|, a set of equations thus obtained has no physical solution. Furthermore, 
the calculations of a triangular Josephson vortex lattice Jl3| , also based on the assumption of the existence of isolated 
Josephson vortices, are at odds with the exact vortex solution of Theodorakis § , valid in the same field range and 
exhibiting full homogeneity in the layering direction. 

As has been shown within the framework of the microscopic theory [^) , the resolution of these contradictions lies 
in the analysis of singular mathematical structure of free-energy functionals of layered superconductors. In partic- 
ular, the system of the Maxwell equations in layered superconductors contains a constraint relation that physically 
constitutes the conservation law for the total intralayer current. According to this constraint relation, the phases of 
the superconducting order parameter (the pair potential) at different layers turn out to be mutually dependent. The 
minimization of the free energy with respect to the phases must necessarily take into account this fact. The neglect 
of mutual dependence of the phases leads to an incomplete set of mean-field equations. In the present paper, we 
elucidate this mathematical issue in full detail. 

Section II of the paper is devoted to the infinite in the layering direction LD model. In subsection II. A, we 
concentrate on exact minimization of the LD free-energy functional. Using general field-theoretical arguments, we 
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prove that the Maxwell equations of the LD model contain an infinite number of unphysical degrees of freedom that 
cannot be eliminated by imposing a gauge condition. We achieve the elimination of these redundant degrees of freedom 
by minimizing the free energy with respect to the phases, taking account of the conservation law for the total intralayer 
current. In this way, we obtain a complete set of necessary and sufficient conditions of an unconditional minimum of the 
LD functional. These conditions constitute a remarkably simple, closed, self-consistent system of mean-field equations 
involving only two variables: the reduced modulus of the pair potential (the same for all the superconducting layers) 
and the phase difference (the same for all the barriers). In addition, we prove that inhomogeneous in the layering 
direction field configurations do not correspond to any stationary points of the free energy. As a result, contrary to 
the prevailing view |||ll],[lO 13 1, the infinite LD model does not admit any solutions in the form of isolated Josephson 
vortices. 

In subsection II. B, we proceed to exact solutions of the mean-field equations of the infinite LD model describing 
major physical effects. We arrive at a new scenario of the flux penetration at H c \. We show that Josephson vortices 
penetrate all the barriers simultaneously and coherently, forming homogeneous field distribution in the layering direc- 
tion (a "vortex plane"). The corresponding lower critical field is H c \ = 2 (nepXj) , where p is the layering period, 

A,/ = (8iTejop) is the Josephson penetration depth, with jo being the critical density of the Josephson current. 
We show that the magnetization exhibits oscillations and jumps due to successive vortex plane penetration. We also 
obtain all well-known limiting cases [the Meissner state, Fraunhofer oscillations of the critical Josephson current, the 
upper critical field H c2 (T), and the vortex solution of Theodorakis] with self-consistent refinements. All these results 
stand in complete agreement with our previous microscopic consideration 

In section III, we consider the finite, both in the layering direction and along the layers, LD model. We show that 
the emergence of additional boundary conditions in this case completely eliminates unphysical degrees of freedom of 
the Maxwell equations and makes minimization with respect to the phases impossible. An explicit solution to the 
Maxwell equations obtained in this section, in contrast to the infinite case, does not preclude the existence of localized 
Josephson vortex configurations. As regards the physical effects, we derive exact analytical expressions for the order 
parameter, the currents and the local magnetic field describing the Meissner state. 

In section IV, we present a brief summary of the obtained results and discuss some theoretical and experimental 
implications. In appendix A, we obtain an explicit solution to the Maxwell equations in the infinite case. We also 
consider a variation of this solution induced by variations of the phases. In Appendix B, we discuss relationship to 
the microscopic theory H . This discussion casts light on the actual domain of validity of the LD model. 



II. THE INFINITE LD MODEL 



In this section, we consider an infinite in the layering direction LD model. One of the dimensions of the system 
along the layers is taken to be finite, although it can be made arbitrarily large. 

We begin by reminding basic features of the LD model Qlfl]. In this model, the temperature T is assumed to be 
close to the "intrinsic" critical temperature T c q of individual layers: 
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< 1. 



(2.1) 



The superconducting (S) layers are assumed to have negligible thickness compared to the "intrinsic" coherence length 
£(T) oc r -1 / 2 , the penetration depth A(T) cx r -1 / 2 , and the layering period p. Taking the layering axis to be x, 
choosing the direction of the external magnetic field H to be z [H = (0, 0, H)], assuming homogeneity along this axis 
and setting h = c = 1, we can write the LD free-energy functional as 
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Here A = (A,,,A y ,0) is the vector potential, continuous at the S-layers: A(np — 0, y) = A(np + 0, y) = A(np,y); 
W z is the length of the system in the z direction (W z — > oo); / n (y) [0 < / n (y) < 1] and <^n(y) are, respectively, the 
reduced modulus and the phase of the pair potential A„(y) in the nth superconducting layer: 

A„(y) = A(T)/„(y)ex P 0„(y), (2.3) 

with A(T) being the "intrinsic" gap [A(T) oc r 1 / 2 ]; H C {T) is the thermodynamic critical field; r(T) = 2a p f l T~ 1 is a 
dimensionless phenomenological parameter of the Josephson interlayer coupling (0 < a p h < 1). The local magnetic 
field h = (0, 0, h) obeys the relation 

h{xy)= 9A^_dA^ 1 y)_ (24) 



A. Exact minimization of the LD functional 

Our task now is to establish a closed, complete, self-consistent system of mean-field equations for the pair potential 



A„ and the local magnetic field h, which is mathematically equivalent to the minimization of (2.2) with respect to 
f n , <j) n , and A. This problem should be approached with a great deal of caution because of singular mathematical 
structure of the functional fl2.2|), resulting from gauge invariance combined with discreteness. Thus, one must take 
account of the fact that variations with respect to (f> n and A are not independent. Moreover, variations with respect 
to <p n at different layers in themselves turn out to be mutually dependent. Unfortunately, these crucial points have 
not been realized in previous literature. To clarify them, we consider partial variational derivatives with respect to 
<p n , and A Xl A y , for mally obtained under the assumption of the independence of these variables. 
As the functional ( |2.2j ) is invariant under the gauge transformation 

<l>n{y) <f>n(y) + 2er?(np,y), Ai(x, y) -> Ai(x, y) + diJ](x, y), i = x,y, 

where rj(x, y) is an arbitrary smooth function of x, y in the whole region (— oo < x < +oo) x [L v x < y < L V 2\, partial 
variational derivatives with respect to </>„, and A Xi A y are related by the fundamental identities 

^ SQld _ d SQ LD | SCl LD 5Vt LD 

8(j> n (y) dy SA y (np,y) 5A x (np + 0,y) SA x (np - 0, y) ' 

Being a consequence of Noether's second theorem, such identities are typical of any gauge theory pifl . They imply 
that the number of independent Euler-Lagrange equations is less than the number of variables. In other words, the 
system of Euler-Lagrange equations contains unphysical degrees of freedo m w hose number is equal to the number of 
Noether's identities. Unusual is, however, an infinite number of identities (|2.5|) . Indeed, in continuum gauge theories 
[such as, e.g., the Ginzburg-Landau (GL) theory of superconductivity] the number of Noether's identities is equal 
to the number of independent parameters of the relevant gauge group. [In the case of superconductivity, we are 
dealing with the electromagnetic one-parameter group U(l).] Thus, by imposing gauge conditions in continuum gauge 
theories, one completely eliminates all unphysical degrees of freedom. By contrast, in the discrete LD theory a single 
available gauge condition cannot eliminate an infinite number of unphysical degrees of freedom resulting from ( |2.5| ). 
The resolution of the problem of the remaining "infinity minus one" unphysical degrees of freedom lies in implicit 
mutual dependence of the variations with respect to the phases cj> n at different S-layers. Below, we demonstrate this 
dependence explicitly. [See relation ( 2.12j ).] 
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To finish with the discussion of ( |2.5| ), we point out that these same identities hold also for the LD model with 
decoupled S-layers [when r(T) = 0]. However, now the num ber of unphysical degrees of freedom is equal to the 



number of physically independent systems [one identity (2.5) per independent S-layer]. A single gauge condition 



completely eliminates the arbitrariness of the Euler-Lagrange equations in this case. 

We start by minimizing with respect to A. Varying ( |2.2| ) with respect to A Xl A y in the regions (n — l)p < x < np 
under the condition 8A x (x, L y i) — SA x (x, L y 2) = yields the Maxwell equations 

dhix vi 

— t^ 1 = ^jn.n-i(y) = 47rj /„(y)/n-i(y) sin (y), (2.6) 

dy 
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dx 
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(2.7) 



where j„,n-i(y) is the density of the Josephson current between the (n — l)th and the nth layers, jo 
r(T)p/16ire( 2 (T)A 2 (T). Minimization with respect to A y (np,y) leads to boundary conditions at the S-layers 
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Equations (p.6|)-(p.q) should be complemented by boundary conditions at the outer interfaces y = L y i,L y2 - As we 
do not consider here externally applied currents in the y direction, the first set of boundary conditions follows from 
the requirement that the intralayer currents vanish at y — L y i,L y 2- 
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(2.9) 



Applied to Eqs. (2_i), these boundary conditions show that the local magnetic field at the outer interfaces is inde- 
pendent of the coordinate x: h(x,L y i) = h(L y i), h{x,L y 2) — h(L y 2). Boundary condit ions imposed on h should be 
compatible with Ampere's law h(L y i) — h(L y \) = AttI obtained by integration of Eqs. (2^) over y, where 
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(2.10) 



is the total Josep hson current. 

Differentiating ( |2.8| ) with respect to y and employing ( |2.6|) , we arrive at the current-continuity laws for the S-layers: 
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These relations may be interpreted as "the Euler-Lagrange equations for the phases" in the s ense that they can be 
formally obtaine d by taking partial variational derivatives with respect to tp n under conditions ( |2.9| ). However, actual 
minimization of ( |2.2j ) with respect to the phases must take account of mutual dependence of 6<p n (y) at different layers, 
as shown in wh at fo llows. [The fact that relations (2.11) follow directly from the Maxwell equations ( |2 . 6[ ) , fl2.8| ) is a 
consequence of ( |2.5| ). Surprisingly, this trivial functional dependence of the current-continuity laws for the S-layers 
on the Maxwell e quat ions has not been pointed out in the previous l iter ature. p5[ ] 

Adding Eqs. ( [2.1l| ), integrating and using boundary conditions ( |2.9| ), we get the conservation law for the total 
intralayer current: 
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This key relation of our consideration has mathematical form of a constraint [ pq | on the derivatives of the phases 
and the y components of the vector potential at different S-layers. Unfortunately, the existence of the constraint 
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relation (2.12) in the system of the Maxwell equations (2.6)-(2.£) has not been noticed in previous publications, hence 
difficulties in establishing a complete set of necessary and sufficient cond itions of an unconditional minimum of ([2.2]) . 
We w ant to emphasize that the fundamental constraint relation (2.12) and its corollaries below [relations ( [2.14 ), 
( 2.15 )1 should not be confused with auxiliary constraint relations imposed on independent variables in the standard 
variational problem of a conditional minimum [T^ ] . All constraints of the LD model appear as a result of singular 
structure of the functional (2J2) itself. (See Refs. flj],[l8| for a thorough discussion of singular field theories of this 

type-) 

According to main principles of the calculus of variations |lq |, to minimize ( |2.2| ) with respect to 4> n , we must first 



eliminate the constraint (2.12). Assuming that f m (y) > 0, where m is an arbitrary layer index, we rewrite (2.12) as 
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Equation (2.13) expresses A v {mp,y) as a function of all 2y ■ ^ should be substituted into (2.2). Now all 54> n (y) 
can be considered as independent. Carrying out the variation under the conditions ([2.S|), we obtain 



/m-i(2/)sin$ mjm _i(?/) - / m+ i(y)sin$ m+ i jm (y) = 0, 



(2.14) 



d_ 

dy 



/»(») 



dcj>n{y) 

dy 



2eA y (np,y) 



d_ 

dy 



fn(y) 



dcj>m{y) 
dy 



2eA y (mp, y) 



r(T) 
2C 2 (T) 



fn(y) [/n-i(2/) sin $ n> „_i (y) - /„+i(y)sin$„ +li „(y)] , n^m. 



Comparing these equations with ( 2.11 ) and integrating with boundary conditions (|2.9|) for n = m yields 
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Since m is an arbitrary layer index, relations (2.14), (2.15) hold for all n = m = 0, ±1,±2, . . . Note that only one 
of the two s ets o f relations (2.14), (2.15) is independent. For example, relations ( |2.14|) can be obtained by inserting 
( |2.15D into ( |2.11| ) and vice versa. In turn, the num ber of independent relations (|2.15|) is exact ly e qual to "infinity 



minus one", because they obey the constraint ( 2.12 ). As expected, the correct minimization of (2.2 ) wi th respect to 
the phase s com pletely resolves the problem of unphysical degrees of freedom contained in Eqs. (2.6)-(2.8). Physically, 
relations ( [2.15|) , which appear already in the case of decoupled layers, minimize the kinetic energy of t he i ntralayer 
currents and, by (2.8), assure the continuity of the local magnetic field at the S-layers. [According to (2.7), h does 
not depen d on x in the barrier regions. Thus, h(x, y) = h(y) in the whole region (— oo < x < +oo) x [L y \ < y < L y 2]-] 
Relations (2.14) constitute stationarity condition s for the Josephson term in (2.2) and assure the continuity of the 
Josephson current at the S-layers as required by ( 2.10| ). 

The above r esult s, in fac t, prove that inhomogeneous in the layering direction field co nfig urations [i.e., those that 
do not satisfy (2.14), (2.15)] do not correspond to any stationary points of the functional (2.2). Consider the variation 
of the solution of (2.(i)-( |2.9| ) for A y in the gauge A x — on an interval (m — l)p < x < mp, induced by a variation 

of the phase at the nth layer. According to (A5), we have SA y (x,y) — /? fa) dS ^ yS> . Such a variation does not 
affect the energy of the magnetic field in ( |2.2| ). If n = m, the variation of the kinetic energy of the intralayer currents 
vanishes, but the first-order variation of the Josephson term is nonzero. If n ^ m, the variation of the Josephson 
term vanishes, but now the first-order variation of the kinetic ene rgy o f t he in tralayer currents is nonzero. These 
first-order variations of ( [2.2j ) vanish if and only if the conditions ( 2.14 ), ( 2.15 ) are fulfilled (i.e., for homogeneous 
field configurations). Unfortunately, this general mathematical consideration unambiguously precludes the existence 
of isolated Josephson vortices ]9|- [Tl]Jl^] in the infinite LD model. It also explains the results of Farid jl2| , who has 
pointed out inconsistencies in a mathematical description of such hypothetical entities. 

It is instructive to look at the incompleteness of the system (2.6)-(2.8) from a slightly different mathematical 
point of view. In the gauge A x = 0, this system reduces to an infinite set of integrodifferential equations ( |A3| ) 
for the phase differences n! „_i (for fixed /„). There are no theorems of existence and uniqueness of a solution to 
an infinite set of such equations. By contrast, for a finite set, describing a finite in the layering direction layered 
superconductor, the existence and uniqueness of a solution can be proved by standard methods of functional analysis. 
The description of a finite layered superconductor implies the specification of boundary conditions on A at the "top" 
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and "bottom" S-layers, whereas the infinite LD model consider ed h ere does not impose any boundary conditions on 
A at x — > ±00. Thus, the arbitrariness contained in Eqs. ( [2.6| )-(2.8) is an intrinsic mathematical property, necessary 
to satisfy additional boundary conditions in the case of the finite LD model. This issue is discussed in more detail in 
section III. 

Minimization with respect to f n is straightforward. Under the condition that df n (L y i), Sf n {Lyz) are arbitrary, we 
get 
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(2.17) 



Equations (|1|-(|J), (|2.16|) and (|2.14j) [or, equivalently, ( |2.15|) 1 (with m — * n), together with boundary conditions 
.9), (2.8) and boundary conditions for h{ y), form a closed, complete set of necessary and sufficient conditions of 
all the stationary points of the functional (2.2). For example, the well-known maximum £Ild = for H = I = 



(the normal state) trivially satisfies these conditions with /„ = 0. The absolute minimum CIld 



Hj{T)V 



(V is the 



volume of the system) for H = I = also satisfies these conditions with <J> n +x,n = and /„ = 1. Complemented by 
the requirement that the Josephson term be a minimum, these conditions become necessary and sufficient conditions 
of all the minima of (2.2) for H ^ 0, I ^ 0, provided that FIld < 0. (For L y 2 — L y i < +00, the Josephson term is 
bounded and thus has both minimum and maximum values.) 

I ndeed , the Josephson term is minimized automatically. T he k inetic energy of the intralayer currents is minimized 
by (2.15). The energy of the magnetic field [the last term in ( |2.2| )] reaches its minimum value for given H and / too. 
This term is non-negative and necessarily has a minimum determined by the condition that its first-order variation 
vanish. (No other stationary points are available.) In the gauge A x = 0, the first-order variation of the magnetic-field 
energy has the form 
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The vanis hing of the volume variation in ( p. 18 ) (the first term on the right-hand side ) is assured by the Ma xwel l 
equations (2.7). The surface variation [the second term on the right-hand side of ( [2.18 )] vanishes by virtue of ( |2.8| ) 
and (2.15). Consider now the condensation energy in ( |2.2[ ) (the sum of the first three phase- an d fiel d-independent 
terms). This energy reaches its absolute minimum for /„ = 1, i.e. when the right-hand side of (2.17) is identically 
equal to zero. The Josephson term and the kinetic energy of the intralayer currents induce spatial dependence and a 
reduction of f n , which increases the condens ation energy. This influenc e is m inimized under the considered conditions: 
the second term on the right-hand side of ( 2.17 ) vanishes according ( 2.1 5| ) and the first term is minimal when the 
Josephson energy is a minimum. 

Thus, we have proved that the above obtained conditions minimize all the terms of the functional (|]^): the 
condensation energy, the Josephson energy, the kinetic energy of the intralayer currents and the magnetic-field energy. 
Any deviation from a solution satisfying these conditions increases all these terms. As a result, the overall LD free 
energy increases, as should be the case for an unc ond it iona l minim um pq|. 

Now we proceed to the simplification of Eqs. (2.6)-(2.8), (2.16) and (2.14) (with m 
field h does not depend on x in the whole region (—00 < x < +00) x [L y ± < y < L y 2\ 
depend on the layer index: 



n). As the local magnetic 
the quantities /„, <J> n! „_i cannot 
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f n {y) = /„_!(y) = f(y), * n +i,n(v) = ^n-i(y) = (2.19) 

The remaining unphysical degree of freedom of Eqs. ( |2.6| )-( |2~8| ), ( |2.14| ), related to the gauge invariance, is eliminated 
by fixing the gauge: 



A x {x,y) = 0, A y (x 1 y) = A(a;,y). 



[Note that dA/dx and d 2 A /dxdy are continuous at the S-layers by virtue of (2 
second set of relations ( 2.19| ) now yields (j> n (y) 



(2.20) 



( ]2~T5| ), and @, (gig).] The 



+ where </>(?/) is the coherent phase difference (the same 
at all the ba rriers), and T](y) is an arbitrary function of y that can be set equal to zer o without any loss of generality. 
From (2.7), employing the continuity conditions for A, dA/dx and relations (2.15), we obtain 
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2ep dy 

Making use of these results, we reduce the functional (|]^) to 
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where W x = L X 2 — L x \. The desired closed, self-consistent set of mean- field equations for the pair potential A n (y) 
and the local magnetic field h(y) takes the form 



A n (y) = Af(y) exp [in<t>{y)\ 
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(2.24) 

(2.25) 

(2.26) 

(2.27) 
(2.28) 

(2.29) 



where h(y) should satisfy appropriate boundary conditions at y — L y i, L y 2 with I = J dyj(y) [see Eq. ( 2.10|) above] 



Remarkably, the c ohere nt phase difference <f> (the same for all the barriers) obeys only one no nline ar second-order 
differential equation (|2.26| ) with only one length scale, the Josephson penetration dept h Xj [Eq. ( 2.27 )1, as in the case 
of the Ferrell-Prange equation for a single junction [Q. [Ma them atically, equation (|2.26| ) is a solvability condition 
for the M axwe ll equations.] Due to the factor / 2 , equation ( [2.26 ) is coupled to nonlinear second-order differential 
equation ( [2.24 ) des cribin g the spatial dependence of the supe rcond ucting order parameter / (the sa me fo r all the S- 
layers). Equations ( |2.25|) constitute boundary conditions for ( 2.24 ). The Maxwell equations ( [2.28| ), ( 2.29| ), combined 
together, yield Eq. ( 2.26| ), as they should by v ir tue o f self-consistency. 

It is important to note that equations (2.23)-( 2.29| ), with an appropriate microscopic identification of r(T) and jo, 
can be consid ered as a lim iting case of the true microscopic equations [Q]. (See Appendix B for more details.) 

Equations ( |2.24[ )-(^.29|), together with ( ]2.22|) , encompass the whole physics of the infinite LD model in parallel 
magnetic fields. They admit exact analytical solutions for all physical situations of interest. These solutions are 
discussed in the next subsection. 
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B. Major physical effects 



1. The Meissner state 



Consider a semi-infinite (in the y direction) LD superconductor with r(T) <C 1, L y \ = 0, L y 2 — > +oo in the external 
fields 



< H < H s = (epXjY 



(2.30) 



In the Meissner state, j(y) — + 0, h(y) — ► for y — * +oo. The requirement that the Josephson term in ( 2.22| ) be a 
minimum means that the density of the Josephson energy should vanish at y — + +oo. This leads to the boundary 
conditions 



^(0)=2epH, ^(+oc) = 0, 0(+oo)=O, /(+oo) = 1. 



(2.31) 



The solution of Eqs. J2.24j ), Q2.26[) , (|2.28|) , ( |2.29| ), subject to ( |2.25| ) and ( |2.3l| ), up to first order in the small parameter 
r(T) has the form 



cf){y) = —4 arctan 



H exp 



(2.32) 



h(y) 



2HH S H s + y/H% - H 2 exp 



H a + y/H% - H 2 + H 2 exp 

j{y) 



2y 



2tt\, 



H a + x //; II- 



(2.33) 



[H s + ^H 2 - H 2 


2 

— H 2 exp 


2y" 
A. ; 


exp 


y 

Aj 




H s + ^H 2 - H 2 


2 r 

+ H 2 exp 




2 



(2.34) 



tf 2 



/(y) = 1 - 4r(T) 



exp 



2j£ 

Aj 



H s + ^/flf - H 2 + H 2 exp 



■2y_ 

A j 



n 2 " 



(2.35) 



The Meissner solution persists up to the field H s = (ep\j) 1 that should be regarded as the superheating field of the 
Meissner state . This fact was established for the LD model by Buzdin and Feinberg [p|. A self-consistent solution of 
the type (p^)-(p5|) was first obtained in the framework of the microscopic theory M. In fields H > H s , only vortex 
solutions are possible. 



2. The lower critical held i/ c ioo- Vortex planes 



Consider now an infinite (in the y direction) LD superconductor with r(T) <C 1, L y \ — oo, L V 2 — > +oo, and 



j(y) -> 0, h{y) -> for y -» ±oo. We are interested in topological solutions of Eqs. Q2.24| ), ( |2.26| ), ( |2.28| ), Q2.29D for 
this situation. The requirement that the Josephson term be a minimum should now be understood as the condition 
that the density of the Josephson energy vanish at y — > ±oo. Thus, the appropriate boundary conditions are 



(-oo) = 0, 



-oo) = ±2tt, 



-f (±oo) = 0, 
dy 



/(±oo) = l. 



(2.36) 
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[Note that aside from </>(+oo) — oo) = ±2ir no other topological boundary conditions are possible. This fact can 
be proved analogously to the well-known case of the sine-Gordon model 0).] 
The desired solutions up to first order in the small parameter r(T) are given by 



4>{y) = ±4arctanexp 



JL 



(2.37) 



h(y) = ±H S cosh " 



JL 

A./ 



(2.38) 



j(y) = T2 jo cosh" 



JL 

A./ 



sinh 



JL 
Aj 



(2.39) 



/(</) = l-4r(T)- 



exp 



2M 
Aj 



1 + exp 



2M 
A., 



These solutions explicitly satisfy the usual conditions of the phase and flux quantization 
rectangular contour T joining the points (— yf>, — oo), (~y Pi +oo) , (+xP> +oo) and 
of the phase along this contour for the "plus" sign in ( |2.37| ) is 



(2.40) 



Indeed, consider a closed 
^■p, — oo). The total change 



A r = dy 



<i0 jv (y) 
dy ~ = 2^iV. 
dy 



Analogously, the total flux through this contour is 



+oo 

*r = Np J dyh(y) = iV$ , 



where $o = is the flux quantum. Thus, the solution with the "plus" sign describes a chain of Josephson vortices 
positioned in the plane y — (one vortex per each barrier). Such a solution was first obtained in the framework of the 
microscopic theory Q and termed "a vortex plane". The solution with the "minus" sign in (2.37) describes a chain 
of Josephson an tivort ices in the plane y — (i.e., "an ant i vortex plane" ). 

By inserting ( 2.37 ) with the "plus" sign and ( 2.40 ) into ( 2.22 ) and comparing the result with the free energy of 
the Meissner state, we derive the lower critical field -ff c ioo, at which the vortex-plane solution becomes energetically 
favorable: 



H 2 H 2 $ ° 

7T 7T TTpAj 



(2.41) 



Note that h(Q) = H s > H cloo . This means that the penetration of Josephson vortices at fields H cloo < H < H s can 
be prevented by a surface barrier, which should result in hysteretic behavior of magnetization [|. Finally, we point 
out that simultaneous Josephson vortex penetration, envisaged by the vortex-plane solution, and hysteresis in the 
magnetization have recently been observed experimentally on artificial low-temperature superconducting superlattices 
Nb/Si ||]. 



3. The vortex state in intermediate fields 



Now we turn to finite-size (in the y direction) LD superconductors with r(T) < 1, -I 



J y2 



W/2, in the field 



range H s < H < H c2oo (H c2 . 
The boundary conditions on 



oo is the upper critical field) and in the absence of externally applied current (1 = 0). 
have the form 



1 # 
2ep dy 



±- 



W 



H. 



(2.42) 
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Under these conditions, the phase difference up to first order in the small parameter H 2 /H 2 is 



<j){y) = 2e P Hy + irN v (H) 



(-1)"- H* 
4 H 2 



[sin (2epHy) — 2epHycos (epWH)] . 



(2.43) 



The constant of integration ttN v (H) accounts for the requirement that the Josephson term in the free energy be a 
minimum. The "topological index" N v corresponds to the number of vortex planes and is a singular function of the 
applied field H: 



N V (H) = 



~epWH~ 






7T 




$0. 



(2.44) 



Here [u] means t he int eger part of u, and $ = pWH is the flux through one barrier. 

By the use of (2.43), we derive the following expressions for the physical quantities up to first order in the small 
parameters r{T) and H 2 /H 2 : 



h(y) = H 



(-!)*' H 2 
4 H 2 



cos (2epHy) — cos (epWH)] 



(2.45) 



j(y) = (-lnh sm(2epHy), 



(2.46) 



"CO 



i - 



(-1) N - cos(2epHy) V2ep((T)H |sin (epWH)\ cosh 
1 + 2 [epC,(T)H] 2 



V2y 



1 + 2 [ep((T)Hf sinh 



w 



V2((T) 



(2.47) 



In the limit W > ((T), \y\ < W/2, equation (2.47) becomes 



r{T) 



{-l) N * cos (2epHy) 
1 + 2 [ep((T)H] 2 



(2- 



The vortex s olutio n (|2.43| ), ( 2.45| ), ( 2.48 ) for N v = 2m (m is an integer) was first obtained by Theodorakis ||. 

From Eq. (2.44) with N V (H) = 1, we derive the lower critical field H ciw in a finite along the layers superconductor 
with W < Aj: 



clW 



(2.49) 



For the magnetization M 



l 



M(iJ) 



/ <^2//i(y) ~ ^ we obtain: 
H 2 [\sm (epWH)\ 



16nH 



1) N " cos (epWH) 



epWH 



(2.50) 



The magnetization (2.50) shows distinctive oscillatory behavior and discontinuities when epWH/it approaches an in- 



teger, i.e when a vortex plane penetrates or leaves the superconductor. For H 3> &o/pW . 



N v ® 
pW 



<H < IN V 



2 J pW 



the LD superconductor exhibits a small paramagnetic effect, i.e. M(H) > 0. (Note that oscillations and jumps of mag- 
netization due to Josephson vortex penetration have been experimentally observed on superconducting superlattices 
Nb/Si d|.) 



4- Fraunhofer oscillations of the critical Josephson current 

Consider the case of a finite-size (along the layers) LD superconductor with r(T) -C 1, — L v \ = L y 2 = W/2 in the 
presence of an externally applied current / in the x direction. The boundary conditions on <p now are 

1 d * f ±™\=H±2.I. (2.51) 



2ep dy \ 2 
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Assuming W <C Aj, we obtain the solution up to first order in the small parameters r(T) and W 2 /X %: 

4>{y) = 2epHy + nN v (H) + tp 



(-l) Nv W 2 _ 2 

(epWH) [sin (2epHy + tp) — 2epHy cos (epWH) cos tp — sin tp] 

4 At 



(2.52) 



T 2 

I(<p,H) = J dyj(y) = -^\sm(epWH)\smtp, 



(2.53) 



h(y) = H 



4 A 2 



1 - 



(epWH) [cos (2epHy + tp) — cos (epWH) cos ip] 



(2.54) 



f(y) = 1 - r(T) 



1 - 



(-1)^" cos(2epHy + tp) 
1 + 2 [epC(T)fl] 2 



V2ep({T)H | sin (epWiT)! cos 95 cosh ^ 



1 + 2 [epCCW 



sinh 



(-l) Ar "\/2epC(r)iZcos(epiyiJ)sin(p sinh 



V2y_ 

C(T) 



1 + 2 [epC(T)fl]- 



cosh 



V2C(T) 



(2.55) 



The phase shift irN v (H) , induced by N v vortex planes, assures the condition of a minimum of the Joscphson energy. 
The field-independent phase shift tp (\tp\ < n/2) parameterizes the total Josephson current I given by ( |2.53 ). Equation 
(2.53) yields the well-k nown Fraunhofer pattern in layered superconductors H§|. Note that the first zero of the 
Fraunhofcr pattern, by ( 2.49 ), corresponds to the lower critical field H cl w- (See Ref. for the explanation of the 
Fraunhofer pattern in terms of the pinning of the vor tex p la nes by t he ed ges of the superconducto r.) In t he ab sence 
of the transport current, i. e., for tp = 0, equations ( 2.52 ), ( 2.54 ), ( 2.55 ) reduce, respectively, to ( 2.43| ), ( |2.45| ) and 



5. The upper critical field H C 2oo{T) 



Here we consider an infinite (in the y direction ) LP superconductor with —L y \ = L y 2 = W/2 — ► +00, subject 
to boundary conditions on the phase of the type (2.42). Supposing that at the upper critical field H = i/ C 2co the 
t rans iti on to the normal phase is of the second-order type, f 2 can be considered as a small parameter, and equations 
(Ol, (E||) become: 



f(y) + C 2 (T) 



d 2 f(y) 
dy 2 



-r{T) [l-cos#y)]/(i/) = 0, 



(2.56) 



d 2 Hv) 

dy 2 



= 0. 



The relevant solution of Eq. (2.57) is 



<j>{y) = 2epHy + irN v {H). 
[Compare with fl2.43Q .]. The substitution of ( |2.58| ) into ( |2.56| ) yields 



(2.57) 



(2.58) 
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(ft fit) 

+ [A(T, H) - (-lf" +1 q(H) cos 2t] f(t) = 0, (2.59) 

lep((T)H] 2 ' 2[e P aT)Hf [e P ^ h H] 2 ' 

where we have introduced a dimcnsionless variable t = epHy and the notation £(0) = £ p h- Hence one gets two 
independent equations: for the odd N v = 2m + 1 (m = 0, 1, 2, . . .) and the even N v = 2m number of vortex planes. 



Both of them have the usual form of the Mathieu equations pi]]. [Note that for N v = 2m Eq. (2.59) is well-known 



The upper critical field H C 2 is now determined by the smallest eigenvalue of ( 2.5S ) 



A{T,H c2oo ) =a (q c ), (2.60) 

where q c = q(H C 2oo), and ao(q) [ao(— q) = ao(q)] is the smallest eigenvalue of the Mathieu equation corresponding to 
the eigenfunctions fN v =2m+i(t) o^ceo(t,q) and /jv„=2m(i) occeo (7r/2 — t, q). [Note that the function ceo(i, q) is strictly 
positive and periodic with the period it.] Explicitly, equation ( |2.60| ) reads: 

~ = 2aph 2 = a Q ( 2* ^ . (2.61) 

[ep£,p h H c2 oo} \[epS, P hH c2 c " 



zoo I 



Equation (2.61) exhibits the well-known 3D-2D crossover of H C 2oc{T) |3|], with the crossover temperature determined 
by r* = 2a p f l . As usual, it is of interest to consider two opposite limiting cases. 
High temperatures, weak fields: r <C 2a p h, -ff C 2oo "C %JOi p hl 'ep£,ph- 
In this 3D regime, 

H c2ao {T) = = 7^^- ( l - ifr) ■ (2-62) 

2y/a ph ept, ph 2^/aph epE, ph \ J-coJ 

The superconductivity of the S-layers is strongly depressed by the vortex planes, which can be seen by comparing 
local maxima / max with local minima / m i n of the order parameter: 

^ = 2V2exp[-2r(T)] < 1. 



Low temperatures, strong fields: 2 (1 — 5^7 J ^ 1) ^c2oo 3> y / &ph/ep£,ph- 
In this regime, 

Hc2oo {T) = ^ ( 1 - — V * = , ^ (2.63) 

This expression diverges for r — * t* — 0. The origin of this well-known unphysical divergence is the unrealistic 
assumption of the LD model of a negligible S-layer thickness. [In the microscopic theory Q], iT C 2ooCO i s finite at any 
temperatures.] The spatial dependence of the order parameter is given by 

(-l) N *r(T) 

f(y) K 1 ~ - r Z7=: TT 72 cos {^pH c2oo y) . 

4 [ep((T)H C 2oo\ 



This spatial dependence is exactly the same as in the case of intermediate fields ( 2.48| ). 



III. THE MEISSNER STATE IN THE FINITE LD MODEL 



Let the LD superconductor occupy the region [L x \ — < x < L X 2\ x [L y \ — < y < L y 2\- The external magnetic 
field H (0 < H < H s ) is again applied along the z axis. The homogeneity along this axis is assumed (W z — > +00). 
The Meissner state realizes under the conditions L X 2 3> A, L y 2 3> A,/, thus it is sufficient to consider the limiting case 
L X 2 — ► +00, L y2 — > +00. 
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This situation is described by the functional (2.2) with a minor change: the summation is now done over n = 
0, 1,2, . . . We assume that r(T) <C 1. Boundary conditions of the type ( ^ ) are supposed to hold at y = 0, and 
h{x 1 0) = H . The presence of an outer boundary at x — is accounted for by the obvious boundary condition 



H -h(+0,y) 



pfi(y) 

2eA 2 (T) 



d(j>o{y) 
dy 



-2eA y {0,y) 



(3.1) 



[Compare with (2.8).] The imposition of the boundary condition (3.1) implies a restriction on variations of A y (x,y): 
they must now satisfy the condition 



5A y {0,y) =0. 

The influence of the boundary at y = must vanish for y — > +oo, hence boundary conditions 



$„+i,„(+oo) = 0, 



-(+oo) = 0. 



For a; — > +oo, we must arrive at the solution of the infinite LD model (2.38)-( 2.40| ), thus 

d</> n (y) 



dy 



2eA y (np, y) — * 0, n — * +oo. 



(3.2) 



(3.3) 



(3.4) 



The minimization with respect to /„ leads to (2.16) and fl2.17|) . Varying with respect to A x under the condition 
SA x (x, 0) = 5A x (x, +oo) = yields the Maxwell equation s (|2.6|) in the regions np < x < (n + l)p (n = 0, 1, 2, . . .). 
Taking variations with respect to A y under the condition ( |3.2[ ), we obtain the Maxwell equations ( |2.7D in the regions 
np < x < (n + l)p (n — 0, 1, 2, . . .) and boundary conditions at the S-layers ( pT§| ) for n = 1, 2, . . . 

Th e general solution of the Maxwell equations, subject to the above formulated boundary conditions, in the gauge 
( 2.20 ), has the form 



A(0,y) = 



r(T) 



4eC 2 (T)/ 2 ( y ) 



dufi(u)f (u) sin^i ;0 (-u) 



1 d<j>o{y) 
2e dy 



A{x,y) 



47rjo J duf n+ i(u)f n (u) sin </)„+!,„(«) + H 
o 



[x - (n + T)p] + 



1 d(j) n+1 {y) 
2e dy 



r(T) 1 
4eC 2 (T)/ 2 +1 (y) 



duf n+ i (u) [/„ (u) sin <f> (u) - f n+2 (u) sin^ n+2 ,n+i(w)] 



(3.5) 



np < x < (n + l)p, n = 0, 1, 2, . . . , 
where the phase differences 0„+i. n obey the solvability conditions 

7/ 

d<pifl{y) 



dy 



87rejoP J dufi(u)fo(u) sin cj) lfi {u) + 2epH 
n 



r(T) 
2( 2 (T) 



fKv) 



dufi(u) [fo(u) sin^i :0 (u) - f2(u) sin 02,1(1*)] 



dufi (u)f (u) sin ^i,o(u) 
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d<j)n+l,n(y) 

dy 



8irej p / duf n+ i(u) f n (u) sin (j) n+1 . n (u) + 2epH 



r{T) 
2C 2 (T) 



duf n+1 (u) [f n (u) sin (u) - f n +i{u) sin</>„ +2: „ + i(u)] 



duf n (u) [f n -i{u) sin 4> n , n -i (u) ~ f„+i(u) sin <f> n+ i tn (u)] 



, 71=1,2, 



(3.6) 



Equations ( p-9 ) assure the continuity of the solution (3.5) at a; = rip (n = 0, 1, 2, . . .). [Compare with Eqs. (A2), (A3) 
of the infinite LD model.] The obtained solution explicitly satisfies the current-conservation law 



87reA 2 ^ 

n=0 



fn(y) 



d<t> n (y) 

dy 



2eA y (np,y) 



duj{u) = 0, 



(3.7) 



where j(y) = lim jn+Xn{y) is the density of the Josephson current given by (2.34). [Compare with the current - 

n — >+oo ' 



conservation law (2.12) of the infinite LD model.] 



Note that, in contrast to the infinite LD model, the minimization with respect to the phases (p n now is not possible. 



Indeed, a variation of the phase at the mth layer, 8<j> m , would induce, by (3.7), a non- vanishing variation of the vector 
potential at x = : 

XA , n .a 1 fm(y) dSKM 

5A ^ V) = Yem^y— 



Such a variation is not allowed by the condition (3.2). This observation explains the role of the arbitrariness contained 
in the system of the Maxwell equations of the infinite LD model, discussed in the previous section. In the infinite 
case, the variations 6A y (±oo,y) were arbitrary, which allowed the system to adjust boundary conditions at x — * ±oo, 
so as to minimize the free energy with respect to the phases. 

We are i ntere sted in the behavior of ( |3.5| ) in the asymptotic region y — > +oo. The second set of the conditions ([3.3]), 
applied to ( |3.6D , in first order in the small parameter r(T) yields 



where 



p \ p H 

In+2,n+l — ( 2 + — I I n +l,n + In.n-1 — = ^' 71 = ^' ^' 



+00 



(3.8) 



I n +l,n=jo / duf n+1 (u)f n (u) sin^„ + i, n (u), n = 0,1,2,... 



is the total Josephson current between the (n + l)th and the nth layers. The local magnetic field inside the barriers 
np < x < (n + X)p (n — 0, 1, 2, . . .) in the asymptotic region y — > +oo is given by 

K+i = h(x, +oo) = 47rJ n+ i,„ + H. (3.9) 

With the help of the quantities h n , equations ( |3.8|) can be rewritten in the form of a recursion relation 



K+2 - ( 2 + ^ I h n+i +K = Q, n = 0, 1, 2, . . . , 



(3.10) 
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subject to boundary conditions 

ha = H, h T 
The solution of ( 3.1C ), ( |3.11 ) is straightforward: 



0, n — > +00. 



H 



P 

2A 2 



P 

4A2 



Assuming p -C A, we get 

h(x, +00) = H exp 



(n+l)p 



np < x < (n + 1) p, 



The intralayer currents are given by 



J„(+oo) 
The order parameters are 



1 

47rp 



(/i„ — h n +i) 



H 



47tA(T) 



exp 



rip 



A(T) 



0,1,2,.. 



0,1,2,. 



/ n (+oo) = 1 - 2e^(WT)i^exp 



2np 



A(T) 



n = 0, 1,2, . 



(3.11) 



(3.12) 



(3.13) 



(3.14) 



(3.15) 



Equations ( 3.13] )- (3.15) describe the Meissner state in t he reg i on [0 < x < +00) x (Aj <C y < +00). In the region 
(A <C x < +00) x [0 < y < +00), the solution is given by ( 2.38 )-( 2.40|) . As in the case of the infinite LD model, the 
upper bound of the existence o f the se solutions is H = H s . Unfortunately, in the region [0 < x < A) x [0 < y < Aj), 
an analytical solution to Eqs. ( 2.16 ) and (3.6) is not possible. 

Equations of the type (2.16), (3.5) and (3.6), subject to topological boundary conditions on tf> n +i. n , in principle, 
describe Josephson vortex configurations. In contrast to the equations of the infinite LD model, these equations do 
not preclude inhomogeneous in the layering direction vortex solutions. 



IV. DISCUSSION 



We have solved the problem of exact minimization of the LD functional in both the infinite and the finite cases. 
We have shown that the LD model belongs to a class of sin gular fi eld theories ]17| , ^8| : the Maxwell equations of 
this model contain constraints [the current conservation laws (|2.12|) , (3.7)1 on the phases and the vector potential at 
different superconducting layers. Such constraints, resulting from gauge invariance combined with inherent discrete- 
ness, are typical of laye re d su perconductors. [See Eq. (26) of the microscopic theory 0].] Unfortunately, the current 
conservation laws (2.12), ( |3.7[ ) were completely overlook ed in previous publications on the LD model. 

By taking into account the current-conservation law ( [2.12 ), we have minimized the free energy of the infinite LD 
mode l with respect to the phases, obtainin g a cl osed, complete, self-consistent system of mean-field equations (2.23)- 
( [2.29 ). We show that relations of the type ( 2.11 ), erroneously regarded as "equations mini mizi n g the free energy with 
respect to the phases" Jl5| |, are, in re ality , mere consequences of the Maxwell equations (2.6)-(2.8). By considering 
non- vanishing first-order variations of (2.2) caused by the variation ( |A5| ) of an inhomogeneous solution to the Maxwell 
equations, we prove that the infini te LD mode l does not admit solutions in the form of isolated Josephson vortices. 

The exact mean-field equations (2.23)-(2.29) contain the whole physics of the infinite LD model in parallel magnetic 
fields. In particular, they reproduce such well-known limiting cases as the Meissner state ||, Fraunhofer oscillations 
of the critical Josephson current the upper critical field H C 2(T) and the vortex state Q in the intermediate 

field regime. All previous calculations of these effects were based on a physical assumption of the homogeneity of the 
solution in the layering direction. We provide a rigorous mathematical justification of this assump tion by proving that 
the homogeneity of the solution is one of the necessary and sufficient conditions of a minimum of (2.2). Moreover, our 
approach allows us to make self-consistent refinements on these results by obtaining exact analytical expressions for 
all physical quantities of interest up to leading order in small parameters, which substantially elucidates the physics. 

We have obtained an exact topological solution (2.37)-(2.4C) to Eqs. d2.23| )-(2.29), describing a chain of Josephson 
vortices (a vortex plane). This solution clearly demonstrates that, contrary to previous suggestions |) 10 1, Josephson 
vortices of the infinite LD model form simultaneously and coherently (one vortex per each barrier) at the lower critical 
field H claol given by (2.41). Successive penetration of the vortex planes at higher fields is accompanied by oscillations 
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and jumps of the magnetization, as described by (2.50). We show that the vortex-plane solutions of the infinite LD 
model persist up to the upper critical field H C 2oo(T) in the whole temperature range. 

Our consideration of t he finite LD model illuminates the role of the boundary effect. Thus, the imposition of 
the boundary condition ( |3.l| ) completely eliminates all unphysical degrees of freedom and makes minimization with 
respect to the phases impossible. The explicit solution (3.5), ( |3.6| ) to the Maxwell equations, in contrast to the infinite 
case, does not preclude the existence of localized Josephson vortex configurations. Making use of this soluti on, w e 
obtai n the first, to our mind, self-consistent description of the Meissner state of the finite LD model [Eqs. ( 3.13Q - 
(E3T51) and ( |2~38[) - (|2~40| )1. 

All the above results stand in full agreement with our previous consideration of layered superconductors j|] based 
on a completely different, microscopic approach. Our discussion (Appendix B) of relationship to the microscopic 
theory Q clarifies a microscopic background of the phenomenological parameters of the LD model and casts light on 
its actual domain of validity. 

As regards the experimental status of the problem, simultaneous Josephson vortex penetration into all the barriers, 
as described in our paper, has been recently observed on artificial low-T c superlattices Nb/Si p2fl . Oscillations and 
jumps of the magnetization, accompanying this penetration, have also been observed Concerning the reported 
observation of localized Josephson vortex configurations in layered high-T c superconductors p3| , the boundary effect 
discussed in section III of our paper may account for this situation. Moreover, the presence of irregularities within the 
layered structure (e.g., stacking faults) can substantially modify the physical picture: such irregularities should serve 
as pinning centers for isolated Josephson vortices. We hope that our exact results will stimulate further theoretical 
and experimental investigation in these directions. 



APPENDIX A: THE EXPLICIT SOLUTION TO THE MAXWELL EQUATIONS OF THE INFINITE LD 

MODEL 



In the gauge A x = 0, the explicit solution of Eqs. (2.6)-(2.8) on the intervals 

(n - l)p < x < np, 7i = 0,±l, n = 0, ±1,±2, . . . , 
subject to boundary conditions (|2.9|) at y = L y \ and h(L v \) = H — 2nl, has the form 



(Al) 



A y {x,y) 



4irj / du f n (u)f n -i(u) sin 4> n , n -i{u) + H -2-kI 

Lyl 



(x — np) + 



l d(j>n{y) 

2e dy 



r(T) 1 
4< 2 (T) my) 



duf n (u) [/„_i(u) sin i(w) - f n+ i(u) sm<j> 

n-f-l,n 

(«)] 



(A2) 



where the phase differences 



>n,n-l — <Pn — <Pn-l 



obey the solvability conditions 



dy 



Snejop J dufn+i{u)f n (u) sin 4> n+ i, n (u) + 2ep(H - 2nl) 



r(T) 
2C 2 (T) 



Pn + M 



(u) - f n +2(u) sin 

(«)] 



my) 



duf n (u) [f n _i(u) sin^> nin _x(u) — f n+1 (u) sin0 

n+l,n 

(«)] 



(A3) 



1G 



This infinite system of integrodifferential equations assures the continuity of the solution (A2) at x — np (n — 
0,±1,±2, . . .). For /„ = 1, equations ( |A3|) reduce to an infinite set of se cond -order non-linear differential equations 
p3| . Unfortunately, an explicit solution for the vector potential of the type ( A2) was not found in previous publications. 

Consider the variation of the solution ( A2) induced by a variation of the phase at one of the S-layers, 8<j) n (y). As there 
is only one constraint on ^ L and the y components of the vector potential at x = np, namely the current-conservation 
law ( 2.12]), s uch a variation affects the solu tion ( A2) only on one of the intervals (|Al[ ), say, (to— 1) < x < mp. Making 
use of ( 2.13| ), we rewrite the solution (A2) on this interval as 



47rj J duf m (u)f„ 

Lyl 



sin0 m)m _i(u) + H — 2ttI 



(x — mp) 



1 d( t>™{y) 1 \- .2, , 

2e dy + Pm{y) ^y) 



1 d(j> n {y) 



2e dy 



A y (np,y) 



(A4) 



In equation (|AJ), all should be considered as independent. Thus, the desired variation is 



5Ay(x,y) = 



i fl{y) dSMv) 

2ef? n (y) dy : 



(A5) 



where n = 0, ±1, ±2, . 



APPENDIX B: RELATIONSHIP TO THE MICROSCOPIC THEORY 

The free energy functional of the microscopic theory [Q, after the minimization with respect to A, has the form 



2 

O = ^p-W x W z J dy 



-f(y) + \f 4 (y) + C 2 (T) 



df(y) 
dy 





\d4>(y)] 




dy 



f(y) 



a£,o 



[l-cos^(y)] f 2 (y) 



+4e 2 ( 2 (T)\ z (T) 



2ep dy 



H 



(Bl) 



where £o is the BCS coherence length, £(^0 aiK1 ^(T) are the GL coherence length and the penetration depth, 
respectively, a is the S-layer thickness, 



3tt 2 
7C(3) 



J dttD{t) < 1, 



(B2) 



with D{t) b eing the tunneling prob ability of the barrier between two s ucce ssive S-layers. The rest of notation is the 
same as in ( 2.23| ). Expression (Bl) applies to the temperature range (2.1) and the S-layer thicknesses meeting the 
condition £ < a < ((T),\(T\ 



Consider the LD limit of (Bl), when a <C p. In this limit, the average kinetic energy of the intralayer currents, 
i.e. the term proportional to a 3 /p 3 , sho uld be dropped. However, the microscopic functional (Bl) does not reduce to 
the corresponding LD functional (2.22) because of the presence of the first order factor a/p. (In the LD model this 
factor is unrealistically taken to be unity.) Nevertheless, as can be easily seen by minimizing ( Bl) wit h respe ct to / 
and (f>, the microscopic mean-field equations in this limit formally coincide with the LD equations (2.23)-(2.29), if one 
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identifies r(T) with the microscopic parameter and the LD quantity j with the microscopic expression for the 

critical Josephson current of a single junction 

Jo = ^|^eiV(0)^A 2 (T), 

where N(0) is the one-spin density of states at the Fermi level. 

The role of the first-order factor a/p becomes evident when one considers the penetration of an external parallel 
magnetic field in the layering direction. As can be shown on the basis of the microscopic equations |2j , the exponential 
falloff of the magnetic field occurs o n the length scale of the effective penetration depth A e // = Ay^T whereas the 



LD model gives A e // = A. [See Eq. (3.13) 
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